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Abstract 

The title equation, where p > 3 is a prime number ^ 7 (mod 8), q is an odd prime 
number and x, y, n are positive integers with x, y relatively prime, is studied. When 
p = 3 (mod 8), we prove (Theorem 2.3) that there are no solutions. For p ^ 3 (mod 8) 
the treatment of the equation turns out to be a difficult task. We focus our attention 
to p = 5, by reason of an article by F. Abu Muriefah, published in this journal, vol. 128 
(2008), 1670-1675. Our main result concerning this special equation is Theorem 1.1, 
whose proof is based on results around the Diophantine equation 5a; 2 — 4 = y n (integer 
solutions), interesting in themselves, which are exposed in Sections 3 and 4. These 
last results are obtained by using tools such as Linear Forms in Two Logarithms and 
Hypergeometric Series. 



1 Introduction 

Diophantine equations of the form px 2 + c — y p , where c is a nonzero integer and p is an 
odd prime, have been studied by several authors. When c = 2™, the case p = 3 was solved 
by Rabinowitz in [27], while Le dealt with the case p > 3 in [24]. The case c = 3™ was 
considered by Abu Muriefah in [2]. Cao [17] treated the cases c = 1 and c = 4" (see also [1], 
[5], [11], [19] for closely related results). We should moreover mention that the equation has 
no solution in positive integers x, y when c = —1, as can be inferred from the work of Nagell 
[26] and Cao [16]. 

The case when c = q 2n , where q is an odd prime, was studied in the recent paper [3] 
of Abu Muriefah. Let us first note that, for fixed n and p > 5, the Diophantine equation 
px 2 + q 2n = y p and, more generally, the Diophantine equation pX 2 + Y 2n = Z p , have at most 
finitely many solutions (x, q, y) and (X, Y, Z), respectively. Indeed, in this case, 5 + 5^ + ^ < 1 
and the claim follows from Theorem 2 of [20]. A main result in the aforementioned paper 
[3], namely, Theorem 3.1, states that the equation 5a; 2 + q 2n — y 5 , has two families of 
solutions given by y = </>3fc,<fofc+i (or ip3k+i, ^3fe+2), k > 1, where 4> k (respectively ip k ) 
is the kth. term of the Fibonacci sequence (respectively the Lucas sequence). However, 
straightforward computations show that the only Fibonacci or Lucas number y < 1000 
satisfying the title equation when p — 5 is y — 21 with x — 410, q — 1801, n = 1, and, 
further, in 5 x 183630 2 + 160201 2 = 181 5 , 160201 is prime and 181 is neither a Fibonacci 
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nor a Lucas number. The same Theorem 3.1 of [3] states that, if p > 3 is a prime ^ 7 
(mod 8) and q is another odd prime, then there are no integer solutions (x,y,n) to the 
equation px 2 + q 2n = y p with (x,y) = 1. The proof of Theorem 3.1 in [3], just before 
its end, contains an obvious, non-rectifiable error, at case 2, when — Wapb 2 is set equal to 
— I6apq 2m although b = ±g J with < j < m. That the said proof is erroneous is also 
pointed out by P.G. Walsh in his review [30] of [3]. One of our aims in the present paper is 
to prove Abu Muriefah's assertion when p = 3 (mod 8); see Theorem 2.3. Our proof, rather 
than rectifying Abu Muriefah's argument (this is probably impossible), goes through totally 
different lines. Unfortunately, our arguments cannot be extended to the case p = 1 (mod 4). 

As we revisited the title equation, we further discovered some new results, like Theorems 
1.1, 3.2 and 4.1 which, we believe, merit one's attention. Moreover, since the powerful 
techniques of sections 4 and 3 are also applicable (after the appropriate modifications) to 
Diophantine equations other than the ones treated in this paper, we thought it useful to 
expose them in some detail, enough for the reader to profit from them. 

As we stated above, one of the main results of this paper is the following 

Theorem 1.1. Let q be an odd prime. If either q ^ 1 (mod 600) or q < 3 • 10 9 , then there 
is no integer solution (x,y,n) to the equation 

5x 2 + q 2n =y 5 , x,y,n>0. (1) 

Otherwise, there exists at most one integer solution (x,y,n) and if it actually exists, then it 

must satisfy the following conditions: 

(i) n < 820 and gcd(n, 2 • 3 • 5 • 7 • 11 • 13) = 1. 

(it) There exists an integer v such that x = 10w(80w 4 — 40v 2 + 1), y = 20v 2 + 1, q n = 
2000v 4 - 200v 2 + l. 

Remark. If the prime q is of the form q = 2000w 4 — 200v 2 + 1 (the first few primes of 
this shape are 1801, 160201, 1245001,4792201, 8179201), then (x,y,n) = (I0v(80v 4 - 40w 2 + 
1), 20v 2 + 1, 1) is a solution to (1) and, according to the theorem, this is the only one (with 
x,y > 0). We have not been able, however, to find a prime q such that the corresponding 
equation (1) has a solution with n > 1. Therefore, we state the following 

Conjecture. // the prime q is not of the form q = 2000?; 4 — 200w 2 + 1, then the equation 
(1) has no solutions. 

The proof of Theorem 1.1 follows from a straightforward combination of Corollary 2.2 
and Theorem 4.1, a second main result of our paper. In turn, the proof of Theorem 4.1 
relies on a third main result, namely, Theorem 3.2 concerning the equation bx 2 — 4 = y n 
which is interesting for its own sake. Indeed, in recent years, important papers are devoted 
to equations of the form x 2 + C = y n . One main strategy for attacking such equations is 
based on the so called modular method which has been successfully applied in quite a number 
of cases; see Chapter 15 (by S. Siksek) in H. Cohen's book [18], the survey article [4] and 
[14] and the references therein. For our equation 5x 2 — 4 = y n , the existence of the trivial 
solution (x, y) = (1, 1) makes the modular approach unsuccessful and prevents us from giving 
the complete solution (x,y,n). Thus, our Theorem 3.2 offers only a partial result which, at 
present, seems to be best possible. 
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2 The Diophantine equation px 2 + q 2n = y p 

The main results of this section are Proposition 2.1 and Theorem 2.3 

Proposition 2.1. Let p > 3 be a prime number ^ 7 (mod 8) and let q be an odd prime 
number. If x, y, n are positive integers with x, y relatively prime such that 

px 2 +q 2n = y p , (x,y) = l,n>0, (2) 

then there exists a rational integer a such that 

(p-i)/2 

,2i + l 



i=0 

and 



± f= E (jm + I U^- 1 ^)^- 1 ^ 2 (3) 



(P"l)/2 



E (S)« p - 2< (-p) (p - 2i - 1)/2 (4) 

Proof. The condition (x, y) = 1 implies p^q and the condition p ^ 7 (mod 8) implies 
that y is odd. 

We work in the imaginary quadratic field K = Q(lo), where ui = \/—p. Equation (2) 
factorizes as (ujx + q n )(—ujx + q n ) = y p and, trivially, the factors in the left-hand side are 
relatively prime. This implies an ideal equation (ujx + q n ) = I p , where I is an integral ideal 
of K. Since the ideal-class number of K is strictly less than p (see page 199 of [21]), the 
above ideal equation implies that I is a principal ideal, therefore we obtain the equation 

„ „ . . \uj if p = 1,5 (mod 8) . , 

w + f-cP, a = aO + b, 0=^ { ^ d8) > (5) 

for some rational integers a and b. In case p = 3 (mod 8) we write the above equation as 
q n — x + 2x8 = (b + a6) p , implying (b + a6) p = or 1 (mod 2). From this, we see that a 
cannot be odd. For, otherwise, we would have 9 P or (1 + 9) p = or 1 (mod 2). But we 
easily check that, for k ^ (mod 3), it is true that 8 k , (1 + 9) k = 6 or 1 + 9 (mod 2), a 
contradiction. Therefore, a in (5) is even and (5) is equivalent to the simpler equation 

xuj + q n = a p , a — aui + b , (6) 

for some rational integers a and b of opposite parity since y = a 2 p + b 2 is odd. Also, it is 
easy to see that (pa, b) = 1. If we put 

= -a, 

then 

_ a p + /3 p q n _a p -/3P 
X ~ 2uj ' T ~ a-P 

and the fact that (a p — f3 p )/(a — f3) is an algebraic integer implies that b divides q n (in Z), 
hence b — ±qi for some j & {0,1, ... , n}. 

At this point we note that the pair (a, —(3) is a Lehmer pair for which 

(a 2 - (3 2 ) 2 = -16pa 2 b 2 . 
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Concerning j appearing in the relation b — ±qi (see a few lines above), we distinguish two 
cases. 

(i) j > 0. Then, in the terminology of [9], (a,—/3) is a p-defective pair. By Theorems 1.4 
and C of [9] it easily follows that p — 5 is the only possibility. Then, by Theorem 1.3 of [9], 
either 20a 2 = (f>k-2e for some k > 3, or 20a 2 = ipk-2e f° r some k ^ 1, where e G {— 1, 1}, 
(0n)«>o denotes the Fibonacci sequence and (ip n )n>o 1S defined by i[) = 2, ipi = 1 and 
V^„ = Vn-i + V , n-2 f or n > 2. It is easily checked that, for every n > 0, ip n ^ (mod 5); 
therefore the second alternative must be excluded. On the other hand, by Theoreme 1.3 of 
[12], a relation of the form <pk = 5z m with m > 1 and z > 1 is impossible, which excludes 
the first alternative as well. 

(ii) j = 0, so that b = ±1. Then, equating rational and irrational parts in (6), we respectively 
obtain the relations (3) and (4). □ 

Corollary 2.2. If the integers x,y,q,n, where (x,y) = 1, q is an odd prime and n > 1, 
satisfy the equation 

5x 2 + q 2n = y 5 (7) 
then [n, 6) = 1, q = 1 (mod 600) and there exists an integer v such that 

x = 10w(80i; 4 - A0v 2 + 1) , y= 20v 2 + 1 , q n = 2000w 4 - 200v 2 + 1 . (8) 

Proof. Applying Proposition 2.1 with p = 5 we obtain ±g™ = 125a 4 — 50a 2 + 1 and 
x = 5a(5a 4 — 10a 2 + 1). Obviously, the minus sign in the first equation is rejected and a is a 
is even. Putting a = 2v in these relations we obtain the first and third relation in (8), and 
then the second relation results immediately. 

We claim that n is odd. Indeed, otherwise (a, q n l 2 ) would be an integral point on the elliptic 
curve defined by Y 2 = 125X 4 — 50X 2 + 1. But this elliptic curve has zero rank and its only 
rational point is (X, Y) = (0, ±1), which forces q = 1, a contradiction. 

We also claim that n is prime to 3. Indeed, let us write the third equation (8) as q n +4 — 5w 2 , 
where w = 20v 2 — 1. If n were divisible by 3, then the last equation could be written as 
(5 2 w) 2 = (5g™/ 3 ) 3 + 500, again forcing q = 1, because it is well known since long (see, for 
example, Table 8 in [25]) that the only integral solutions (X,Y) to Y 2 = X 3 + 500 are 
(X,F) = (5,±25). 

Finally, we show that q = 1 (mod 600). First, we write the third equation (8) as q n — 
200w 2 (f 2 — 1) + 1, which shows that q n = 1 (mod 600). Let r be the order of q modulo 600. 
Then r divides n, and since y(600) = 160 and n is odd, we obtain r = 1 or r = 5. The 
latter case cannot hold, for, otherwise, 5x 2 = y b — (g 2 ™/ 5 ) 5 , which has no proper solutions by 
Theoreme 2(2) of [22]. We therefore conclude that r = 1 and g = 1 (mod 600), as claimed. 

□ 

Theorem 2.3. Letp,q be odd primes withp = 3 (mod 8) andp > 3. Then, the Diophantine 
equation 

px 2 + q 2n = y p (9) 
has no positive integer solutions (x,y,n) with (x,y) = 1. 

Proof. By the relation (3) we see that the equation (9) implies 

q 2 

Let us consider the polynomial 

,, , _ (cox + 1)P - (lux - 1)P 
j[x) .— 
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Clearly, this is a polynomial in Z[x] of degree p — 1, with leading coefficient — p(p +1 )/ 2 and 
constant term 1. 



First Claim: The polynomial /(a;) factorizes over Q[x] into two relatively prime 
polynomials fi(x), $2{x) € Z[x], each of degree (p — l)/2. 

Proof: Let £ be a primitive p-th root of unity, i.e. a root of the p-th cyclotomic polynomial 
$(x) 



c p 1 + • • • + x + 1 = (x p — l)/(x— 1). Let also 5 be a primitive root mod p. Observe 
that the field L = Q(C) contains ui. Indeed, it is a well-known fact that the Gauss sum 
G = YTtZi (p) C* satisfies G 2 = (— l)( p_1 )/ 2 p. Therefore, we can assume u = G. Below we 
will use the well-known fact that, in the field L we have the factorization into prime ideals 
(p) = (A) p_1 , where A = 1 — For (3 £ L we will write w(/3) to denote V\(f3), the exponent 
of A in the prime factorization of /3; and for b e Q we will write v(b) to denote v p (6), the 
exponent of p in the prime factorization of b. 
We have 

p-i 

f(x) = U{(u J x + i)-( k (u J x-l)}, 

k=l 

from which it follows that the roots of f(x) are exactly the following: 

i + C fe 



U) 

P 



1 -c k 



k = l,...,p-l. 



Therefore, 



Let us put now 



fix) 



-V 



(p+1)/2 



p-1 

n 

fe=i 



X 

V 



i-C A 



/l(^) 



/ 2 (x) = 




_ p (p+l)/4 



so that /(x) = /i(x)/ 2 (x). We now show that the polynomials /,(x) have rational coefficients. 
The Galois group of the extension Q(C)/Q is cyclic generated by the automorphism <r, 
defined by a(() = ( 9 . Since uj e Q(C) \ Q> we must have <t(uj) 7^ w, therefore, ct(w) = -co. 
Consequently, for the typical root of /i(x) we have 



i + C fl 
i-C s 



to 

p 



1 + C 9 



i-C 



,2j + l 



1 _ 



(10) 



where 



——+.7 (mod——) 



and we choose 




if <j<H^3 

if^<i<^ ' 
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so that f runs (exactly once) through all values 0, 1, ... , as j runs through these values. 
Consequently, the coefficients of the polynomial fi(x) are fixed by a, which implies that they 
belong to Q; and similarly for f2(x). Actually, the coefficients of fi(x), f 2 (x) are integers 
and we prove this as follows. 

First, we show that the absolute value of the constant coefficient of both f\(x) and f2(x) is 1. 
Indeed, let bi be the constant coefficient of fi(x). We already know that bi € Q. Moreover, 
multiplying the right equalities (10) for j = 0, . . . , (p — 3)/2 and then the resulting products 
in the two sides by — p(p +1 )/ 4 , we obtain b 2 = b\. But, b\b 2 is equal to the constant term of 
f(x), which is 1. Therefore 1 = b\b 2 = b\, from which bi = b 2 = ±1. 
Let us put now 

g(x)=x"- 1 f(-), g 1 (x)=x^- 1 » 2 f \(-) , g 2 (x) = x^ 1 )/ 2 ^-) , 

XXX 

i.e. these are the reciprocal polynomials of fix), fi(x) and f 2 (x), respectively. Since f(x) — 
fi(x)f 2 (x), we also have g(x) = gi(x)g 2 (x). Since the constant term of f(x) is 1, g(x) has 
leading coefficient 1; and since f(x) has integer coefficients, so does g(x). Therefore, the roots 
of g(x) are algebraic integers. Analogously, the polynomials gi(x) have leading coefficients 
equal to ±1, their coefficients are rational numbers and their roots, being roots of g(x), are 
algebraic integers. Therefore, these polynomials have coefficients in Z; consequently, the 
same is true for the polynomials fi(x), as claimed. 

Now, observe that fi(x) and f 2 (x) have no common roots, therefore they are relatively prime. 
Second Claim: Let Res(/i,/ 2 ) be the resultant of the polynomials /i,/ 2 - Then 
Res(/ 1 ,/ 2 )=±2^ 1 ) 2 /V p2 - 1)/8 . (11) 
Proof: We use the symbol Disc to denote the discriminant. We have 

Disc(/) = Disc(/i/ 2 ) = Disc(.f 1 )Disc(/ 2 )Res(.f 1 ,/ 2 ) 2 . 

By the right-most equality in (10) and the comments following it we see that f 2 (x) = iff 
fi(—x) — 0, hence 

Disc(/) = Disc(/i/ 2 ) = Disc(.f 1 ) 2 Rcs(/ 1 , / 2 ) 2 . (12) 

Calculation of Disc(,f ): 

'uj 1 + w 1 + C j 



Disc(/)= P ( 2 p- 4 )(f +i )/ 2 n 

l<i<j<p-l 



p 1-0 p 1-0 



2 



1 \ (p-2)(p-1)/2 / i r i r j\ 

=P (p+i)(i.-2) _±\ n " 



„(p-l)(p-2) (p-2)(p+3)/2 TT ( C'-C 1 
l<Kj<p-l vv s A 



(i-0)(i-o) 



l<i<j<p-l 

2 



The right-most product in the last equality is a unit times A ( p 2 \ therefore, w(Disc(/)) = 
{p- l)(p- 2)(p + 3)/2 - (p- l)(p- 2) = (p- 2)(p 2 - l)/2. But, since Disc(/) is a rational 
integer, it follows that 

Disc(/) = ±2 (p-i)(f- 2 yp- 2 )(P+i)/ 2 . (13) 
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Calculation of Disc(/i): 

( uj l + C 9 ' M w i + c 9 



v\ 2 



Disc(/i)=^- 3 ^ +i )/ 4 n 

0<i<j<(p-3)/2 



v i-C 9 p l-C 



23 



1 \ (p-3)( P -l)/8 / ~ (/ . g ™ 2 3 

:p (H-l)(p-3)/4/Zi) JJ I rJl ' 



p / A1 , V (i-C 9 )(i-C 9 ') 

^ 7 0<i<j<(p-3)/2 \ v s /v s 

= _ 2 ( p - i )( p - 3 )/v p_3)(p+3)/8 n ' ;,< ' ^ 



0<i<j<(p-3)/2 



(i- C ^)(i-C 9 ^) 



The right-most product in the last equality is a unit times A ( p 3 )/ 4 , therefore, w(Disc(/i)) = 
(p-l)(p-3)(p + 3)/8-(p-l)(p-3)/4 = (p-l)(p-3)(p+l)/8. Since Disc(A) is a rational 
integer, it follows that 

Disc(/i) = ±2(f- 1 )(f- 3 )/ 4 p(f- 3 )(P+ 1 )/ 8 . (14) 
Now the relations (12), (13) and (14) imply the validity of the relation (11). 

Third Claim: Among the integers /i(a),/ 2 (a) one is equal to ±1 and the other 
is equal to ±q n . 

Proof: By Bezout's identity, there exist polynomials hi(x) 1 h 2 {x) e Z[x] (both of degree 
< (p - l)/2) such that 



/n^)/!^) + ^w/^x) = Disc(/i,/ 2 ) - ±2(f- i ) 2 /yf 2 



-l)/8 



We make the substitution x <— a in this equality. By /i(o)/2(a) = /(a) = ±g n and the fact 
that (g, 2p) = 1 (cf. beginning of the proof of Proposition 2.1), it follows that exactly one 
fi(a) is equal to ±1 and the other is equal to ±q n . 

Fourth Claim: If fi(a) — ±1 for some i e {1,2}, then a = 0. 

Proof: Let us put fi(x) = c + C\x + • • • + c r x r , where r = (p — 1) /2. We already know that 
c = ±1 and c, € Z for all i. By the very definition of the polynomial /i, its roots are 

1 + C j 

where S* is a complete set of quadratic residues mod p. 

Since $(—1) = 1 the numerator 1 + 7 is a um t and it follows easily that 

w&) = -(p+l)/2. 

Thus, for fc > 0, 

p+lp-l-2fc (p-l)(p+l) fc(p+l) 



w 



(<*)>- » ; + 



2 2 4 2 

Then, for k > 1, we have v(cfe) = w(cfc)/(p — 1) > and therefore 

p|c fc) fc = l,...,(p-l)/2. (15) 
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Moreover, for k > 2, we have v(cfc) > (p + l)/(p — 1) > 1, hence v(cfc) > 2 and we can write 
therefore 

fi(x) = c + c\x (mod p 2 i 2 ) . (16) 

Next, we prove that 

v(ci) = l. (17) 

First, note that /^(a;) = fi(—x), which results from the fact that the polynomial f\ is of odd 
degree, and the polynomials f\ and f 2 have opposite leading coefficients and roots (cf. just 
before the relation (12)). These observations imply that f 2 (x) = c — c\x + c 2 x 2 — c 3 x 3 + ■ ■ ■ . 

y.2 W 



Another observation is that c\ 7^ 0. Indeed, since f(x) = fi(x)f 2 {x), the coefficient of x" m 
f(x) is 2c c 2 

non-zero. 

A third fact - which is a bit more than an observation is that 



is — p 



c\. On the other hand, by the initial definition of f(x), the coefficient of x 2 
which is odd, because p = 3 (mod 4). Therefore, C\ is odd; in particular, it is 



(18) 



If we prove this, then, in combination with the relation (15) and the fact that ci 7^ 0, we 
will conclude that v(ci) = 1. 

Proof of (18): Let gi(x) be, as before, the reciprocal of the polynomial fi(x). Then c\ 
is equal, up to sign, with the sum of the roots of gi(x). But the roots of gi(x) are the 
reciprocals of the roots of fi(x), i.e. they are equal to £J 
(remember that S is a complete set of residues mod p) , 



|ci|< Vp±^ 

jes 



Since 



7r(p-2j) 
2p 



1 + 

< -|, it follows that 



jes 



tan 



TTJ 



i = l,...,(p- l)/2. Therefore 



1 



tan 



*-(p-2j) 
2p 



tan 



7r(p-2j) 
2p I 



> ^b-2j|, hence, 
1 



Note that, as j runs through the set S, the numbers |p — 2j\ are distinct modp, for, if 
\p - 2ji\ = \p - 2j 2 \ (mod p) with j 1} j 2 € S and ji ^ j 2 , then, necessarily, j 2 = -ji, 
which implies that —1 is a quadratic residue modp, a contradiction. Therefore, the set 
{\p— 2j\ : j G S} is a subset of {1, . . . ,p— 1} with cardinality (p— l)/2. It is clear, therefore, 
that 

(P v^ /2 l 3 , p. 



E 



1 



P-2j| 



< 



1 



fe=i 



from which we obtain 



2pyp /3 

7T V2 



log 



p-1 



This upper bound for |ci| clearly implies |ci| < p 2 , as claimed. 

Final step of the proof of Theorem 2.3: By our third claim above, /i(a) or /2(a) must be 
±1. Since /2(a) = — a), we may suppose that /i(a) = ±1, i.e. c + cia + c 2 a 2 + - • • = ±1. 



Remember that cq = ±1, therefore, cq + C\a + C2" + 



±Cq. The — sign implies 



±2 + cia + c 2 a 2 



0, clearly impossible, in view of (15). The + sign implies 



Ci<z + c 2 a 2 + • • • . If a 7^ 0, then, taking also into account (16), we obtain = c\ (mod p 2 ) 
which contradicts (17). This forces a = and then, by (3), q n = ±2, a contradiction. □ 
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3 The equation 5a: 2 — 4 = y 



The third relation (8), written as q n — 5(20w 2 — l) 2 — 4, naturally leads to the study of the 
more general equation 

5x 2 -4 = y n , \x\,y>l, n>2 (19) 

in the integer unknowns x, y, n, where x and y have not, of course, the same meaning as the 
x, y in equation (7). 

First, let n be even. It is well-known that the positive integer solutions of 5X 2 — 4 = Y 2 
are given by X = F 2 k+i, Y = L 2 k+i for k > 0, where F denotes the Fibonacci and L the 
Lucas sequence; notice that k = gives y = 1 which is excluded. Since it is known that the 
only Lucas number which is a pure power is L% = 4 ([13], Theorem 2), it follows that the 
only solution (x, y, n) of the equation (19) with even n is (2, 2, 4). 

From now on we suppose that n is odd. 
If n = 3, then (25x) 2 — 500 = (5y) 3 . It is well known that the only integral solutions (X, Y) 
to Y 2 = X 3 + 500 are (X,Y) = (5, ±25), corresponding to y = 1, which has been excluded 
by hypothesis. Hence we may assume that gcd(n, 3) = 1, in particular n > 5. 
If x is even then 5x 2 — 4 = —4, 16, 12 (mod 32) implying n < 4, which has already been 
excluded. Hence x and y are odd. 

Now we work in the field K = Q(#), where 9 = y/E. From now on and until the end of 

the paper we view K as embedded into the real numbers with 9 ^ V5 = 2.2360679 The 

ring of integers in K is I = {(x + y9)/2 ; x, y G Z with x = y (mod 2)}, e = (1 + Vb)/2 is 
the fundamental unit. In K unique factorization holds. Throughout this section, for a G K , 
a' will always denote the algebraic conjugate of a. We factorize the equation (19) over the 
field K 

Hx 2 -4= (x9 -2)(x9 + 2). (20) 

If p is a (rational) prime divisor of y, then 5x 2 — 4 = y n implies that p is odd and, clearly, 5 
is a quadratic residue mod p. It follows that p splits in K and p = ±1 (mod 5). Therefore y 
factorizes in I as y = irir' , where we can choose n > (then ir' is also positive). Notice also 
that y n = 1 (mod 5), hence y = 1 (mod 10) (remember that y and n are odd) and y > 11. 

Without loss of generality, we assume that x is positive. Since x is odd, x6> + 2 and x6 — 2 
are coprime with x6> > j/™/ 2 > ll 5 / 2 . Hence, there exists fc G Z such that x(9 + 2 = e k Tr n . 
Writing fc = fc + fc! with -(n-l)/2 < fci < (n-l)/2, we have x6» + 2 = ^(eV)" = e fel 7if , 
where 7Ti = e e ir. The conjugate relation is — x6 + 2 = s'^n'™ and summing the two relations 
we get 

£ fcl 7r™ +e" £ V 1 " =4. (21) 

We have 7Ti = u + v9 or 7Ti = (u + v9) /2, where u, v G Z are unknown and in the second case 
uv is odd. Then, for fixed n and fci we obtain from (21) 

T kl («, w) := e fcl (« + «0) n + e'* 1 (u - v9) n = 4 or 2™+ 2 , (22) 

where, in the second case, uv is odd. Note that the left-hand side of (22) is a homogeneous 
polynomial in 7,[u,v] of degree n, hence the relation (22) implies a Thue equation. Since 
T_kA u i v ) = (~ l) fel ^*;i ( u 7 — w )i it suffices to consider the Thue equations T^u^v) = ±4 
and Tfe 1 (u, u) = ±2"+ 2 with fci = 0,1,..., (n — l)/2, where, in the second equation, uv 
is odd. Moreover, since the degree of the form is odd, we can ignore the minus sign 
in the right-hand sides. Using the above Thue equations we will prove that there arc no 
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solutions (x, y, n) to (19) with n € {5, 7, 11, 13}. Actually, we will show that for these values 
of n the Thue equations T^ 1 (u, v) = 2 n+2 with uv odd, and (u, v) = 4 are impossible 
for all fci = 0, 1, . . . , (n — l)/2. For every n as above, the method is practically the same. 
However, as one can guess, the case n = 13 is somewhat more complicated; so we briefly 
expose this case in order to illustrate how we work. Numerous Thue equations of degree 
n arise. A practical method for the solution of such equations has been developed since 
long by Tzanakis and de Weger [28] which later was improved by Bilu and Hanrot [8] and 
implemented in Pari (http : / /par i. math. u — bordeaux.fr/ and Magma [10], [15]. We use 
either of these packages to solve the Thue equations that arise. 

We assume now that n = 13 and we consider all fci's in {0, 1, ... , 6}. 
fci = 0: Since T n (u,v) is reducible, our equations arc treated by elementary means; no 
solutions arise. 

fci = 1: Both equations T\(u, v) = 4, 2 15 are easily solved. 

fci = 2: The congruences T 2 (u,v) = 4, 2 15 (mod 13 2 ) are impossible. 

fci = 3: The equation T^iu, v) — 2 15 with uv odd implies solvability of the congruence 
T 3 (x, 1) = (mod 2 14 ). But, as it is easily checked, this congruence has no solutions. The 
equation T 3 (u, v) = 4 remains. Since T 3 (u,v) = 4u 13 + we multiply by 2 11 and we 

obtain a Thue equation w' 13 + 65u' 12 v H h 320000000u 13 = 2 13 , whose only solution is 

(u 1 , v) — (2, 0) which we obviously reject. 

fci = 4: Now, T 4 (u,v) = 7u 13 H h 234375v 13 . On multiplying by 7 12 we obtain monic 

Thue equations with right-hand sides 4 • 7 12 and 2 15 7 12 . No solutions are returned. 

fci = 5: Similarly to the case k\ = 2, both congruences T 5 (u 7 v) = 4, 2 15 (mod 13 2 ) arc 

impossible. 

fci = 6: All coefficients of T 6 (u, v) are even and \Tq{u, v) = 9u 13 +260u 12 vH h312500w 13 = 

Tg(u, v), say. We thus have the Thue equations Tg(w, v) = 2 14 with uv odd, and Tg(u, v) = 2 2 . 
The first equation implies solvability of the congruence T$(x, 1) = (mod 2 13 ) with x odd, 
which is impossible. For the second equation we are obliged to multiply by 3 24 in order 
to obtain a monic Thue equation, as required by both Pari and Magma. The resulting 
equation is treated with some "effort" by Magma and no solutions are returned. On the 
other hand, Pari after several hours was still "struggling" , so we gave up. 

The computational difficulties arising above, when fci = 6 show the limitation of the 
method and, indeed, for n = 17 the computational difficulties for the solution of the resulting 
Thue equations, at present, seem to be insurmountable. 

Summing up our results so far, we have the following theorem. 

Proposition 3.1. There are no solutions (x,y,n) to the equation (19) with n divisible by 
at least one of the primes 3,5,7, 11 or 13. The only solution (x,y 7 n) to the equation (19) 
with even n is (±2,2,4). 

Computing a first upper bound for n. We now fix a solution (x, y, n) of the equation 
(19), where, in view of Proposition 3.1, we can assume that n > 17. Obviously, we can also 
assume that n is prime. Based on the few observations just after the equation (20), but 
relaxing the condition x > 0, we see that there exists a set V consisting of (unordered) sets 
{7r,7r'} such that 7r > 0, 7T7t' = y and, if {^1,^1} and {tt2 7 tt' 2 } are distinct elements of V, 
then 7T2,7T2 are non- associated to both 7Ti,7Ti. 

We modify V as follows: Let {ir,ir'} € V . There exists precisely an m € Z such that 
e m < \/£y/7r < e m+1 . The last relation is equivalent to e 2m - 1 ir 2 < y < £ 2m+1 ir 2 . On 
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putting e m TT = 7Ti we obtain 

-^-L < y/y < TTiy/e , or, equivalently, J— < m < ^fey. (23) 
Note that tx' x = (— l) m e~ m ir', so that 7Ti |7r^ | = y. On multiplying the first relation (23) by 

y y 

\ir[\ we get — F < |7ri | < vV^i hence W - < |7rJ| < y/ye. The last relation combined with 

V e V £ 

the second relation (23) implies maxl^i/l^l , Itt^. I / 1" 71 "! I } < e anci , certainly the left-hand 
side of the last inequality is > 1. We make the substitution tt •<— tti or tt <— \ir[\ according 
as 7Ti / j 7r j_ | is > 1 or < 1, respectively. In this way, an "adjusted" set V\ replaces the set V 
containing elements {7r,7r'} such that, 

tt>0, 7r|7r'|=y, 1 < ir/\ir'\ < e . (24) 

Now, in view of the relation (20) and the fact that the two factors in the left-hand side are 
relatively prime, we must have an ideal equation (2 + x9) = (ir) n or (2 + x9) = (tt')™ for some 
{7r,7r'} <G Vi, and then (—2 + x9) = (n') n or (—2 + x9) = (n) n , respectively. By choosing the 
appropriate sign for x we may assume that (2 + x9) = {n) n , from which it follows that 

2 + x0 = as k ir n for some k e Z and a e {-1, 1}, (25) 

and 

x Q-2 = ^f— = ?^^=ae- k W\ n . (26) 
x9 + 2 <7S k ir n 11 v ' 

By (25) and (26) we obtain 

WM" 1 ge **" 1 x6 + 2 1 4 (27) 
\\tt'\J ae- k \Tr'\ n x9-2 x9 - 2 ' ( 1 

We have hx 2 = y n + 4, from which \x\9 > y n / 2 . 
Now we put 

k'l 

A = 2fcloge-nlog — , (28) 



7T 



so that A = \og(e 2k (^JT^j ) and now ' b y ( 27 )> 



A 4 4 4_ 4.0001 

|e ~ 1 ~ \x9 - 2\ ~ \x\9 - 2 < y™/2 _ 2 < ~^72~ ' 

Notice that the right most side is less than 5.63 • 10 -9 in view of the fact that y > 11 and 
n > 17. Therefore, 

On the other hand, since the ideals (ir) and (it') are distinct, ir/\n'\ is not a unit and, 
consequently, A ^ 0. Thus, 

0<IM<^ <») 

and 

Tl 

log I A| < --logy + 1.3963. (30) 
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Now we compare k and n that appear in the linear form A. Since log e and log ^ are both 
positive and, by (29), |A| is very small, it follows that k is negative. By (29), |A| < 7.6 x 10~ 6 , 
therefore, in view also of (24), 

|*| = -k = + -!L_ log JL < 7.9 x 10" 6 + lo ge = 7.9 x 1(T 6 + J , 

21oge 2 loge | vr' | 2 loge 2 

hence 

l*l<2" ( 31 ) 

Next, we consider the algebraic number r\ :— appearing in A. This number is a root of 
the polynomial 

(ttX - |vt'|)(|vt'|X - tt) = yX 2 - (tt 2 + tt' 2 )X + y = yX 2 - (a 2 ± 2y)X + y e Z , 
where a = it + n' eZ. From this we easily see that 

Hv) < \ (logy + loge). (32) 

Finally, we are ready to calculate a first upper bound for n using Corollary 2 of [23] . In view 
of the relations (32) and (24) it is easy to estimate the quantities that are involved in that 
corollary. Choosing the parameter m that appears in the corollary equal to 20, and taking 
into account that max(2|fc|, \n\} = n (cf.31)), we easily find that, if n > 15100, then 

log | A| > 78.8 flog n + log lo S^ + lo S £ + 1 + 38 \ ( log + log e) 
V logj/ + loge ) 

This, combined with (30), gives 

(loff v -\- los e 4- 1 \ tl 

logn + log , , +0.38 (logy + loge) - -logy + 1.3963 > 0. (33) 

logy + loge J ' 2 

Since y > 11, we easily check that the inequality (33) can hold only if 

n<2.2xl0 4 . (34) 



Proving that solutions with "small" y cannot exist. Now we go back to our equation (19) 
and we assume that (x, y) is a positive solution. It is easily checked that this positiveness 
restriction does not prevent us from obtaining again the relations (27) and (29). We write 
the last inequality in the following shape: 



log?7 



loge 



2k 
n 



< 



4.0402 



n log e • y"/ 2 ' 



r] :-- 



The right-hand side is, obviously, less than l/(2n 2 ), which shows that 2k/n is a convergent 
to the continued fraction expansion of log n/ loge and, moreover, the denominator of this 
convergent is less than 10 5 , in view of (34). Let do, Oi, (X2, • • ■ be the partial quotients and 
Po/qo,Pi/qi,P2/q2, ■ ■ ■ the convergents to that expansion. Let h be the first subscript such 
that qh > 10 5 . Then, 2k/n = p m /q m for some m £ {0, . . . , h — 1}. We have now 



< 



log?? _ Th 
loge q, t 
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hence, 



1 



(Oi+l 

from which it follows that 



+ 2)r 



< 



log rj 



loge 



2k 
n 



< 



4.0402 



n\og£ ■ y n / 2 



4.0402(A + 2)n>loge-y 



n/2 



A := max{a ,ai, . . . ,a h } . 



(35) 



For every y = 1 (mod 10) with y < 3 • 10 9 and for every it as above (there are 2 m such 77' 
s, where m is the number of rational prime divisors of y), we compute 77 and the continued 
fraction expansion of the real number log 77/ loge, and we check the validity of the relation 
(35). These computations can be performed with the routines of either Pari or Magma. 
We stress the fact that an ordinary precision is sufficient since the denominators of the 
checked convergents have at most 10 decimal digits. The whole task took around 30 hours 
of computations with Pari in a usual PC; with Magma it would take more time. It turns 
out that, except possibly if n < 11, this relation is not satisfied. But we already know that 
n > 17, hence we conclude: 



No solutions (x, y, n) to (19) exist with n > 17 and y < 3 • 10 9 . 

Obtaining a smaller upper bound for n. Now, we know that y > 3 • 10 9 (this is very 
important!) and we apply Theorem 2 of [23] to our linear form A. In the notation of that 
theorem, we choose p = 15.7 and /1 = 0.57. We obtain a new lower bound for |A| which, as 
before, we combine with (30) to obtain a complicated relation in which the only parameters 
that are present are y and n. Since y > 3 • 10 9 , we check that n must necessarily be at most 
6404. Falling from the upper bound (34) to n < 6404 is already a considerable improvement, 
but we can do much better. We turn to the main theorem of [23], namely Theorem 1. Its 
application is somewhat complicated, as one has to choose appropriately various parameters 
p,LL,Ri,R,2,Si,S2,K,L, but it's worth the trouble! The strategy is the following: Once we 
have the lower bound y > 3 • 10 9 and an upper bound for n, we choose our parameters above 
in order to obtain a smaller upper bound for n. We repeat the process with this reduced 
upper bound and new parameters, and so on, as in the table below. In the table we omit 
the values of K and L because we always choose 



L = RiS\ , K = max{ 



R2S2 — 1 
i?l«Si 



,2}. 



The choice of the remaining parameters requires experiments; we cannot expose a systematic 
strategy for choosing them. 

Table 1: The choice of parameters in Theorem 1 of [23] 



p 


n 


Ri 


Si 


R2 


S 2 


starting ub for n 


reduced ub for n 


11 


0.6 


1 


20 


400 


110 


6404 


5802 


11 


0.6 


1 


20 


300 


100 


5802 


3956 


11 


0.6 


1 


20 


211 


80 


3956 


2226 


10 


0.57 


1 


15 


205 


79 


2226 


1949 


9.95 


0.57 


1 


15 


198 


75 


1949 


1783 


10 


0.555 


2 


8 


195 


75 


1783 


1713 


10 


0.56 


2 


8 


186 


66 


1713 


1452 



continued to the next page 
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Table 1 (continued from previous page) 



p 


P 


Ri 


Si 


i?2 


s 2 


starting ub for n 


reduced ub for n 


10 


0.555 


2 


7 


184 


63 


1452 


1358 


10 


0.555 


3 


5 


178 


60 


1358 


1251 


10 


0.555 


4 


4 


173 


60 


1251 


1216 


11.01 


0.6 


4 


4 


172 


49 


1216 


1113 


10.09 


0.59 


4 


3 


170 


48 


1113 


1054 


11 


0.6 


6 


2 


165 


46 


1054 


1003 


10.9 


0.6 


7 


2 


162 


46 


1003 


980 


15.7 


0.57 


4 


3 


149 


42 


980 


901 


15.555 


0.565 


4 


3 


149 


42 


901 


890 


15.3 


0.565 


4 


3 


149 


42 


890 


885 


16.05 


0.58 


4 


3 


145 


41 


885 


878 


16.05 


0.58 


4 


3 


145 


41 


885 


878 


16.045 


0.579 


4 


3 


145 


41 


878 


877 


16.045 


0.579 


6 


2 


142 


41 


877 


859 


16.05 


0.578 


6 


2 


142 


41 


859 


858 


16.051 


0.577 


6 


2 


142 


41 


858 


856 


16.052 


0.576 


6 


2 


142 


41 


856 


855 


16.053 


0.575 


6 


2 


142 


41 


855 


853 


16.054 


0.5744 


6 


2 


142 


41 


853 


852 


16.054 


0.5743 


6 


2 


145 


40 


852 


849 


16.042 


0.5788 


6 


2 


144 


40 


849 


848 


16.042 


0.5788 


6 


2 


160 


35 


848 


825 


16.03 


0.575 


6 


2 


160 


35 


825 


820 



No further reduction of the upper bound of n was possible. 



Summing up all the results of Section 3 we obtain the following theorem. 

Theorem 3.2. Any integer solution (x,y,n) of the equation 5x 2 — 4 = y n with y > 1 and 
n > 2 must satisfy the following: 

(i) n < 820. 

(ii) The prime divisors of n are > 17 and < 811. 
(in) y > 3 • 10 9 . 

Remark. In recent years, the so called "modular approach" to certain types of Dio- 
phantinc equations -the Fermat equation being one of them- turned out to be very succesful; 
see, for example, S. Siksek's "The modular approach to Diophantine equations", Chapter 15 
in [18]. Our equation 5x 2 — 4 = y n resembles the Lebesgue-Nagell equation x 2 + D = y n , to 
which the modular method applies succesfully in most cases; see [14]. However, as mentioned 
in [14], the method is not succesful when D = —a 2 ± 1 because, in that case, there exists an 
obvious solution valid for every n. In the case of our equation we face a similar situation: 
the existence of the solution (x, y) = (1, 1) for every n makes the application of the modular 
method "hopeless", according to S. Siksek (private communication). 
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4 The equation 5a: 2 — 4 = y n when y is prime 



The main result of this section is the following 

Theorem 4.1. Let q be an odd prime. Then, for the solutions (x,n) of the equation 

hx 2 -A = q n , x>0,n>0,n^2 (36) 

the following are true: 

(i) If q ^ 1 (mod 10), no solutions exist. 

(ii) If q < 3 • 10 9 , no solutions with n > 2 exist. 

(Hi) If (q +4)/5 = □, then (x,n) — (yj 1) is the only solution. 

(iv) If (q + 4)/5 7^ □, then at most one solution exists. 

(v) No solutions exist with n > 820. 

(vi) No solutions exist with n divisible by a prime from the set {2, 3, 5, 7, 11, 13}. 

The proof of this theorem follows from a straightforward combination of Theorem 3.2, 
already proved in Section 3 and Proposition 4.2, below. Therefore, the present section is 
essentially devoted to the proof of this proposition. 

As noted in the beginning of Section 3, the third relation (8), written as 5(20v 2 — l) 2 — 4 = 
q n , led us to the more general equation (19) for which Theorem 3.2 holds. In this theorem, y 
is general and not necessarily prime as the equation 5(20v 2 — l) 2 — 4 = q n would suggest. In 
this section, however, we will add the extra restriction that the unknown y in the equation 
(19) be a prime, say y — q, and we will prove the following theorem. 

Proposition 4.2. If q is an odd prime, then the equation 

5.x 2 = q n + 4 , x, n positive integers, n odd (37) 

has at most one solution if (q r + 4)/5 is not a perfect square and exactly one solution, namely, 
(x, n) — (y/ (q + 4)/5, 1) if (q + 4)/5 is a perfect square. 

Remark: It is easy to see that the relation (37) implies q = 1 (mod 10). 

Proof. The proof of Proposition 4.2 will be completed in three steps. 
Step 1: The gap between two solutions of (37). This step consists in proving that, if two 
solutions (x,n), and (x',n') exist, with n' > n, then n' must be "very large" compared to n; 
see (42). We need first the following result. 

Lemma 4.3. Let x be a positive integer and assume that 

2 + xV5 = C , (38) 

where £ is an algebraic integer in Q(V5) and a is an integer > 1. Then, £ = 1+ 2 V ^ , a — 
3,x = l. 

Proof of the lemma. There are two possibilities for £: (I) Either £ = tt^/E with b, c odd 
integers, or (II) £ = b + ca/5 with 6, c arbitrary integers. 
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After expansion of the right-hand side of (38) we obtain 



b a + 5c 2 



ua-2 



+ ... = 



2 a+1 in case (I) 
2 in case (II) 



It follows from this that, if a is even, then an odd power of 2 is a square mod 5 which is 
impossible. Therefore, a is odd and 



b a + 5c 2 



Ld-2 



-I 1- 5— c a 



a 

a- 1 



2 a+1 in case (I) 
2 in case (II) 



Case (II) is impossible. Indeed, note that in the left-hand side all exponents of b are odd and 
all exponents of c are even, hence b > 0. Also, b divides 2, hence 6=1 or 2. If 6=1 then 
an obviously impossible congruence mod 5 results; and if 6 = 2 then 2 a < 2 which implies 
a = 1, contrary to the hypothesis. 

In case (I) we have, as before, 6 > 0, 6 is odd and 6|2 a+1 . Hence, 6=1 and we have 



2 a+1 = 1 + 5 



c 2 + 



+ 5" 



a 

a-1 



c a - x > 2 



a+1 



where the last inequality is strict for every c, if a > 5 and for every c with |c| > 1 when 
a = 3. Thus, to avoid the contradiction we must conclude that a — 3 and |c| = 1 from which 
it easily follows that c = 1 and £ = (1 + \/5)/2. This completes the proof of Lemma 4.3. 

Weput6> = (l + v / 5)/2, 0' = (l-V5)/2. These are the roots of the polynomial x 2 — x — 1 
and is the fundamental unit of the ring of integers of Q(0). In general, for any a € Q(9) 
we denote by a' the conjugate of a under the isomorphism 9^9'. 

Assume now that (x, n) is a solution to equation (37). Then (2 + a;\/5)(2 — xV5) = — q n 
and it is clear that the factors in the left-hand side are relatively prime as algebraic integers 
of Q(6>). Also, every (rational) prime dividing q factors into two distinct prime ideals. It 
follows then that there exists an algebraic integer a with norm ±q such that the following 
ideal relation is true: (2 + xv / 5) = (c)™- Then, for some reZwe have the element equation 
2 + x^fh = ±9 r a n and since we can assume without loss of generality that a > 0, we finally 
get 

2 + xVb = 9 r a n along with the conjugate relation 2 — x\fh = 6' r a' n . 
Combining the last two relations we obtain 



0<5:= ^ - 



9 r a r 



4 1 

WE < 2 



Then, i < 1 - 6 = -(^) r (^)" < 1 and in view of the inequality | log(l - x)\ < \x\(l + \x\ 
(valid for \x\ < 1/2) we obtain 



-r log 



nlog 



<*(! + <*), 6 



< 



2 + xV5 2 + Vg" + 4 q n / 2 ' 



(39) 



Now, let (x',n') another solution to (37) with n' > n, n' odd and x' > (hence, x' > x). 
Exactly as before we have a relation 2 + x'\fh — 9 r a n for a convenient r'eZ and 



-r' log 



n' log 



<5'(1 + S'), 6' 



< 



2 + Vr+4 q 



n'/2 



< 5. 



(40) 



1G 



Putting u = \og\0/0'\ = log((3 + \/E)/2) and eliminating the term log \a'/a\ from the in- 
equalities (39) and (40) we get | - rn' + r'n\u < n'5{\ + 5) + n5'(l + 5') < 2n'6(l + 5), 
i.e. 

2n' 

\-rn'+r'n\< — S(l + 5) . (41) 
u 

The left-hand side in (41) is non-zero. Indeed, in the opposite case we would have 7^ = 77 = 
(say) ^- with (ri,n\) = 1. Then, r = ar\,n = an\, r' — br\,n' — bn\ for some positive odd 
integers a, b with a < b and, moreover, 2 + x'\/5 = (8 ri a ni ) b . By Lemma 4.3 we conclude 
that x' = 1, contrary to the fact that x' > x > 1. 

We conclude therefore that the left-hand side of (41) is > 1, from which it follows that 

„'>Hj-i(i + *)-!, (42) 

which shows that, the larger solution n' is "far away" from the smaller solution n; specifically, 
it is of the size of q n ^ 2 . This fact will play an important role below. 

Step 2: Application of Hypergeometric Polynomials. At this second step we adapt to our 
equation the method of F. Beukers in [6] and [7]. As a result we prove Lemma 4.4 below, after 
which the final step for the proof of Proposition 4.2 is not difficult. In that method one uses 
as a tool the hypergeometric polynomials, the properties of which we remind immediately 
below. 

Given the real numbers a, (3, 7 where 7 is not zero or a negative integer, we define the 
hypergeometric function (with parameters a, /?, 7) 

F ia , A 7 , „ _ 1 + "A z + f; i ^ + i)---(° + t-iWH-i)-- ; tf + t -i) ^ 

7 ^fc! 7(7 + l)---(7 + fc-l) 

which converges for every complex number z with \z\ < 1 and, in case that 7 > a + j3, it 
also converges for z = 1. Let n 2 > n\ > be integers. Put n = m + n 2 and define 

G{z) = F(-n 2 - 1/2, -m, -n, z) , H(z) = F(-m + 1/2, -n 2 , -n, z) . 

By the definition of G it is easy to see that 

^-K-^^) (;)"'(-.)', 

which, in particular, shows that, for any real number z < 0, G(z) is positive. We will use 
the following properties: 

1. G(z) and H(z) are polynomials in z of degrees ni and n 2 , respectively. Moreover, the 
polynomials (^)G(4z) and ( r ™)_ff(4z) have integer coefficients. 

2. \G(z) - (1 - z) 1 /2 jff ( z )| < G(l)|z|™ +1 for \z\ < 1. 

3. G(l) < G(z) < G(0) < 1 for < z < 1. 

4. If G*(z) is the polynomial resulting from G(z) when ni,n 2 are respectively replaced 
by ni + 1, n 2 + 1 and H*(z) is defined analogously, then 

G*(z)H(z) - G(z)H*(z) = cz n+1 

for some non-zero constant c. 
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"2 + 1 



for any z. 



For the proof of the first four properties see Lemmas 1,2,3 and 4 in [6]. For the proof 
of the fifth property see relation (1.10), page 226 of [29]. Now we arc in a position to prove 
the main result of this step. 

Lemma 4.4. Let (x,n) be a solution to (37), where, as always, x > and n > 1 is odd; we 
assume, moreover, that q n > 600. Let r,s be positive integers such that q n > 2 6+4s / r and 
me the positive real number v by means of the relation 



q nv = 2.007 x (4.03g n ) r / s . 
Finally, let N — q n ' where n' > n and let y be any integer. Then, 



AT 1 / 2 



1 



0.27 J (1+u)/2 

^ n n(3+v/2) q 



Proof of the lemma. Let n 2 > n\ be positive integers which will be specified later and 
m = ni + n 2 . Put z = —q~ n . Then, \Az\ < 1 so that G(Az) and H(Az) are meaningful. By 
properties 2 and 3 of the polynomials G and H we have 



G(iz)-H(iz){l+-?l* 



m+1 



' q n 



m+1 



hence 



m )G(4z)-( m )-< ^ 



G(Az) - 



nij 



H(Az)- 



,n/2 



< 



ft) 



n\) \q n 



m+1 



By property 1 and the fact that G(x) > for any negative real number x, 



m )G{Az) 
mj q" 



for some positive A e Z 



and similarly, 



Then, (43) implies 



H(Az) = - 
nij q 



B 



Bx y/5 



Q 



nn 2 (jn/2 



< 



for some B £ Z. 



m+1 

— I , from which 
q" 



BxV5 



Aqn(n 2 -ni+l/2) 



< 2 r 



gTini 2 2 (m+l) 2^m+l 



Now, let us put e 



jvT72 



so that, from the above inequality we have 
Bx 



jyl/2 Aqn(n 2 -m + l/2) 



1 / 2 3m+1 

< 7! V + A(z«(«2+i) j ■ 



Let A = . Then, 



1/2 



< 



(43) 



(44) 



(45) 
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Now comes the moment to choose ni,n 2 - First we choose n\ to satisfy 

1 



r , r , 2s 
-A < m < -A H 

s s s 



(46) 



We must keep in mind that there are exactly two consecutive positive integers in the interval 
[rX/s , (rX + 2s — 1) /s]; this is a simple exercise. Choose now n 2 by setting n 2 = n\ + A > m 
and remember that m = n 1 + n 2 = 2m + A. Moreover, we will need below that the left-hand 
side of (44) be non-zero. In the next lines we show that we can choose n\ in such a way that 
this requirement be satisfied. 

Suppose that for the smaller integer n\ in the interval [rX/s , rX + 2s — l)/s] the left-hand 
side of (44) is zero. Then, we can repeat the above process with := n\ + 1 in place of 
n\ (n'j still belongs to this interval), n' 2 := n 2 + 1 in place of n 2 and m' := n' x + n' 2 in place 
of to, so that the polynomials G and H will be replaced by G* and H* respectively, and 
the integers A, B by some other integers, say, A*,B*. Then, we will obtain an inequality 
analogous to (44), namely, 



y 



B* 



N l/2 A*q n< - n 2- n 'i+ 1 / 2 ) 



1 



23m +1 



< V5 I 6 + A*q n ( n 2+ 1 ) 



If the left-hand side were again zero, then we would have B/A = B* /A* (note that n 2 — n' x = 
n 2 ^ n i), which would easily imply that z = —4/q n is a zero of the function G* ■ H — G ■ H* 
and this contradicts property 4 of the polynomials G, H. We conclude therefore that for at 
least one integer m satisfying (46), the left-hand side of (44) is non-zero and from now on 
we assume that we have selected such anm. 

We now rewrite the term J 4gr™( Tl2_ni + 1 / 2 ) appearing in the left-hand side of (44). We 
first observe that (45) implies q n I" 1 < q"( A+1 / 2 ) which shows that g™( 2A + 1 ) = q n q 2 ^ for some 
non- negative integer fi. Consequently, on putting q^ = A n (a positive integer), we have 

Aq n(n 2 -n 1+ l/2) = Aq n(\+l/2) = Aq n' /2 Aq = AqAN 1/2 

Going back to (44) , we get 



23m+l 



— (e 



> 



> 



y 



Bx 



jyl/2 Aq n(n 2 - ni + l/2) 
1 1 



y 



Bx 



N 1 / 2 AoAN 1 / 2 



\A Ay-Bx\ 
AqIAIN 1 / 2 



^o|^4|^V 1/2 |yl| g «( A +i/2) 



from which 



e\A\q ni - x+1/2) + 2 3m + 1 q -"(™i+ 1 /2) > ^ _ 



(47) 



We estimate separately the second summand in the left-hand side of (47). By the hypothesis 
on the lower bound of q n and (45) we have 

o3m+l o6ni+3A+l o6ni+3A+l 

f ^ f ^ _f _ 23A+l-4sni/r <- 23A+1-4A _ 2 1 ^ A < 1 

q nni q nri\ — 2(6+4s/r)ni — ~ — ' 

which shows that the second summand in the left-hand side of (47) is < q~ n l 2 < 60CT 1 / 2 . 
This shows that the first summand in the left-hand side of (47) is larger than 5 1 / 2 — 60CT 1 / 2 > 



19 



2.195. Then, remembering also how A has been defined and using property 5 of the polyno- 
mial G, we get 



2.195 < e\A\q n( - x+1 ^ = e^+A+i/2) (™ \ \G(-A/q n )\ 



2 \ "2 + 1 



< 



n(m+A+l/2) . 2 m-l A , _M = eQ n/2 n(m+\) . \ ( <± + — 



< |q"/ 2 g "("i+ A ) x 2.007™ (since g" > 600) 

= i q n/2 q nX(l+m/X) x 2 QQ 7 A(l+2ni/A) 

< £ n/2 nA(f + 25=l+l) x 2 ^(^ + ^+1) 

- 2 y y 

= i (? ™/ 2 g ™( 2s - 1 )/ ;s x 2.007( 4s ~ 2 ) /s • (q n ( 1+r / s ) x 2.007 1+2r//s ) A 

< i g n /2 g T.(2«-i)/« x 2 .007 (4s ~ 2)/s • (2.007 • (4.03g") r/ V l ) A 

= £g f n/2 g n(2«-i)/a x 2.007 (4s ~ 2)/s • q"( 1+ ^ A (by the definition of i/) 

< |g n /2 g n(2«-l)/a x 2 . 07 4 . q n(l + v)\ _ 

from which we immediately get 

0.27g" ( -3 + ^ < eg " A ( 1+ ^ < e{Nq n ) iv+1 ^ 2 

(the right-most inequality being true because q nX < {Nq 71 ) 1 / 2 in view of (45)), and hence 

the claimed lower bound for e = | '^{f 2 — 1|. This completes the proof of Lemma 4.4. 

Step 3: Completion of the proof of Theorem 4-2. Assume that, if (q + 4)/5 is not a perfect 
square there exists a solution to equation (37) and if (q + 4)/5 is a perfect square there 
exists a solution to this equation besides the obvious one which results from the relation 
5(y / (q + 4)/5) 2 = q + 4. Thus, in both cases, our assumptions in particular imply that there 
exists a solution (x , no) with n > 1, hence, by Theorem 3.2, we must have q > 3 • 10 9 . Let 
(x, n) be the least solution to equation (37). In order to prove the theorem, we will assume 
that a larger solution (x',n') to (37) exists and we will arrive at a contradiction. 
We put N = q n , so that N + 4 = 5x' 2 , from which we get 

4 ( \/5V ^ ( VEx' , ( VoV ( VFR \ ( V$x' ,\ 

J. ] j> 2i 



N [N 1 / 2 J yN 1 / 2 J [N 1 / 2 I \ y/N J [N 1 / 2 
therefore 

„ V$x' 2 

We apply Lemma 4.4 with y = x', r = 1, s = 2; then it is easy to check that v < 0.7178 and 
by the conclusion of the lemma and the last displayed inequality we get 

27V- 1 > - 1 > 0.27 X q -n{5+u)/2 N -{l+u)/2 ^ 
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licncc 

(1-^) / , 5 + v , 5.627 + 1/ , 
— - — m log < log(7.408) H — n log g < n log g , 

from which 

, 5.627 + u 
n < — n . 

1 — v 

On the other hand, recalling that u = log((3 + V5)/2) and 5 = 2+ ^„ +4 , we have in view of 
(42), 



»<>^-.(i +j) - = H.a±#±3ii + 



2 v ; 2 4 V 2 + V?" + 4 

= u _ (2 + V^+4) 2 u _ (2 + g"/ 2 ) 2 12 x (2 + g"/ 2 ) 2 
8 6 + V?" + 4 8 6 + g™/ 2 ' 6 + g"/ 2 

Combining this lower bound for n' with (4), we get the following relation: 

0.12x( 2 + g "7< 5 - 627 + V 
6 + q n / 2 1 - v 

By the definition of v, {q n ) u = 2.007 x (4.039™) 1 / 2 . Solving for v and substituting into the 
above inequality we obtain 

Q27nlo gg + 7 12 (2 + 9 "/J 7 = log(2 . 007 .4.030- 5 ). (48) 
0.5nlogq-7 6 + q n / 2 ' 6V ; V ; 

However, in view of the large size of q we easily check that (48) is not satisfied and this 
contradiction proves that the solution (x',n') cannot exist. □ 
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